Search for: All records

Abstract We present genome assemblies for 18 snake species representing 18 families (Serpentes: Caenophidia): Acrochordus granulatus, Aparallactus werneri, Boaedon fuliginosus, Calamaria suluensis, Cerberus rynchops, Grayia smithii, Imantodes cenchoa, Mimophis mahfalensis, Oxyrhabdium leporinum, Pareas carinatus, Psammodynastes pulverulentus, Pseudoxenodon macrops, Pseudoxyrhopus heterurus, Sibynophis collaris, Stegonotus admiraltiensis, Toxicocalamus goodenoughensis, Trimeresurus albolabris, and Tropidonophis doriae. From these new genome assemblies, we extracted thousands of loci commonly used in systematic and phylogenomic studies on snakes, including target-capture datasets composed of ultraconserved elements (UCEs) and anchored hybrid enriched loci (AHEs), as well as traditional Sanger loci. Phylogenies inferred from the two target-capture loci datasets were identical with each other and strongly congruent with previously published snake phylogenies. To show the additional utility of these non-model genomes for investigative evolutionary research, we mined the genome assemblies of two New Guinea island endemics in our dataset (S. admiraltiensis and T. doriae) for the ATP1a3 gene, a thoroughly researched indicator of resistance to toad toxin ingestion by squamates. We find that both these snakes possess the genotype for toad toxin resistance despite their endemism to New Guinea, a region absent of any toads until the human-mediated introduction of Cane Toads in the 1930s. These species possess identical substitutions that suggest the same bufotoxin resistance as their Australian congenerics (Stegonotus australis and Tropidonophis mairii) which forage on invasive Cane Toads. Herein, we show the utility of short-read high-coverage genomes, as well as improving the deficit of available squamate genomes with associated voucher specimens. 

Abstract Divergence dating analyses in systematics provide a framework to develop and test biogeographic hypotheses regarding speciation. However, as molecular datasets grow from multilocus to genomic, sample sizes decrease due to computational burdens, and the testing of fine-scale biogeographic hypotheses becomes difficult. In this study, we use coalescent demographic models to investigate the diversification of poorly known rice paddy snakes from Southeast Asia (Homalopsidae:Hypsiscopus), which have conflicting dates of origin based on previous studies. We use coalescent modeling to test the hypothesis thatHypsiscopusdiversified 2.5 mya during the Khorat Plateau uplift in Thailand. Additionally, we use ecological niche analyses to identify potential differences in the niche space of the two most widely distributed species in the past and present. Our results suggestHypsiscopusdiversified ~ 2.4 mya, supporting that the Khorat Plateau may have initiated the diversification of rice paddy snakes. We also find significant niche differentiation and shifts between species ofHypsiscopus, indicating that environmental differences may have sustained differentiation of this genus after the Khorat Plateau uplift. Our study expands on the diversification history of snakes in Southeast Asia, and highlights how results from smaller multilocus datasets can be useful in developing and testing biogeographic hypotheses alongside genomic datasets.